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Abstract 

We characterize numerically the regime of anticipated synchronization in the cou- 
pled FitzHugh-Nagumo model for neurons. We consider two neurons, coupled uni- 
directionally (in a master-slave configuration), subject to the same random external 
forcing and with a recurrent inhibitory delayed connection in the slave neuron. We 
show that the scheme leads to anticipated synchronization, a regime in which the 
slave neuron fires the same train of pulses as the master neuron, but earlier in time. 
We characterize the synchronization in the parameter space (coupling strength, 
anticipation time) and introduce several quantities to measure the degree of syn- 
chronization. 
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Synchronization phenomena is a fascinating subject that has attracted a lot 
of attention in the last years [1]. Recently, Voss[2,3,4,5] has shown that it is 
possible to synchronize autonomous dynamical systems in a master-slave con- 
figuration in such a way that the slave system can actually anticipate (i.e. 
predict) the trajectory of the master system. This result is surprising at first 
sight for two facts: the dynamics of the master is not modified by the presence 
of the slave, and the slave integrates its equations of motion at the same speed 
than the master does. This remarkable phenomenon is achieved by the intro- 
duction of appropriate delay lines in the dynamics of the slave system. More 
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Fig. 1. Schematic diagram of two model neurons coupled in a unidirectional config- 
uration, subjected to the same external stimulus and with a feedback loop (with a 
delay time r) in the slave neuron. 

precisely, one of the schemes devised by Voss considers master x(t) and slave 
y{t) (vector) dynamical systems, whose dynamics follow the general form: 

*(*) = f(x(t)) 

y(i)=f(y(t))+K[x(t)-y(f-r)]. 



The function f (x) defines the dynamical system under consideration, K is the 
coupling strength matrix and r is the delay time in the feedback loop of the 
slave system. As stated before, these equations admit the "anticipated" mani- 
fold y(t) = x(£ + r) as a (structurally) stable solution [2]. This has been shown 
to be possible even in systems in which the dynamics, being chaotic, is highly 
unpredictable. Implementations of this result have been demonstrated theo- 
retically [6] and experimentally [7] in unidirectionally coupled laser systems 
as well as in electronic circuits [5]. Some understanding of the anticipation 
mechanism can be achieved by the study of simple maps [8,9]. 

We have recently extended this result [10] by considering non-autonomous 
systems in which the dynamics is subjected to the effect of an external per- 
turbation. Namely, we consider dynamical equations as: 

±(t) = f (x(t)) + I(t) 

y(t)=f(y(*))+I(0 + K[x(t)-y(t-r)]. 



where I(t) is an external input acting on both the master and slave systems 
(see scheme in Fig. 1). Remarkably, although the manifold y(t) = x(t + r) 
is no longer an exact solution of the equations (2) (except in the case of a 
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periodic input I(t) = I(t + r)), it will be shown that the slave can actually 
predict some interesting part of the dynamics of the master. 



We have considered as our dynamical system x a simple model for neuronal 
response. Neurons are classical prototypes of excitable systems: their response 
to an external perturbation is highly non-linear and depends on its magnitude 
and timing. If the perturbation is small the system evolves back to the steady 
state; but if the perturbation exceeds a certain threshold, the system fires 
a pulse-like spike (action potential). Following the onset of the excitation, 
there is an interval during which another perturbation does not induce a new 
pulse (refractory period). Real neurons are complicated non-linear systems 
involving a large number of variables. Nevertheless, the essential features of 
their excitable behavior can be captured with a much-reduced description. 
The FitzHugh-Nagumo model provides the simplest representation of excitable 
firing dynamics and it has been widely used as a prototypic model[ll,12]. 

In this paper we study numerically the anticipated synchronization of two 
identical FitzHugh-Nagumo neurons, unidirectionally coupled, in the presence 
of a common external random forcing (see the schematic diagram shown in 
Fig. 1). The model equations are: 

x\ = ~xi(xi - a)(x! - 1) - x 2 + I(t) 

I 3 ) 

x 2 = - bx 2 ) 



Vi = - a)(Vi - 1) - 2/2 + /(*) + K[ Xl (t) - Vl (t - r)] ^ 

2/2 = e(yi - by 2 ) 



where x=(xi,x 2 ) are the variables associated to the master neuron, y= (2/1, 2/2) 
are the variables associated to the slave neuron, and a, b, and e are constant 
parameters, n controls the strength of the coupling and r is the delay time 
associated to the feedback loop in the slave neuron. Note that only the fast 
variables x±, y\ are coupled. The external common forcing lit) is a Gaus- 
sian random process of mean I and delta-correlated in time (white noise): 
([/(*) - I ][I(t') -Io]) = DS(t -If), D being the noise intensity. 

By choosing the mean value of the noise Jo just below the threshold of the 
excitable system, a highly complex dynamics is observed. Spikes develop at 
random times in a completely unpredictable manner, see Fig. 2. The same 
figure shows that although the exact details of the master dynamics are not 
reproduced by the slave, still it manages to anticipate the response of the 
master by firing its pulses just before the master does. 

We now proceed to quantify the degree of synchronization. To this end we have 
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Fig. 2. Trains of spikes obtained from numerical simulations of Eqs. (3-4). The 
parameters are a = 0.139, b = 2.54, e = 0.008, I = 0.03, k = 0.1, r = 4, 
D = 2 Ah x 10~ 5 . Left panel: Spikes of the master x\(t) and slave yi(t) neurons. 
Notice that the slave neuron makes an error around t ~ 1000 (4th pulse) in firing 
when the master does not. Right panel: detail of an anticipated spike. The solid line 
is the pulse of the master and the dotted line is the pulse of the slave. 

measured the relative number of errors, R, made by the slave when anticipating 
the dynamics of the master. Notice in Fig. 2 that the slave neuron occasionally 
fires an extra spike, which does not correspond to a spike fired by the master, 
but every spike fired by the master has a corresponding anticipated spike fired 
by the slave. Thus, an error is defined as a pulse in the slave that has no 
corresponding pulse in the master. 




Fig. 3. Relative number of errors R in the parameter space (k, r). The white region 
represents a region where the number of errors is larger than R = 0.1. 



Figure 3 displays R in a gray scale in the parameter space (k, t). The dark 
(white) region represents a region where good (bad) synchronization occurs. 
In order not to miss too much detail, errors larger than R = 0.1 have been 
uniformly plotted as white, while black indicates R = 0, and the gray levels 
run between these two values. Two different synchronization mechanisms are 
present in Fig. 2. The first one appears for very low coupling intensity (the 
black region near the ordinate axis k — 0). This is not a regime of anticipated 
synchronization, but it corresponds to the synchronization of trajectories by 
common random forcing [13] which leads simply to x(t) = y(t). 

Beyond this regime of synchronization by common random forcing, a finite 
value of the coupling k is required in order to achieve anticipated synchro- 
nization. However, a very large value of the coupling worsens the quality of 
the synchronization. This counter-intuitive result can be explained since the 
dynamics of the slave system becomes chaotic for such large values of the 
coupling and cannot follow the dynamics of the master. The existence of min- 
imum and maximum values for the coupling in order to exhibit good antici- 
pated synchronization agrees with what was previously found in autonomous 
chaotic systems [2] and in linear maps [9]. 

To quantify the anticipation time, we have computed the mean value (t) and 
standard deviation a of the time difference t™ — tf, where t™ are the times 
when the master neuron fires a pulse, and tf are the times when the slave 
neuron fires the corresponding pulse (hence the erroneous pulses fired by the 
slave are not taken into account). The data shown in the next figures are the 
result of averaging over a few thousand spike events. 

Figure 4(a) plots the mean anticipation time (t) as a function of r for different 
values of the coupling k. The results for large k fall mainly on the line (t) = r 
corresponding to the anticipated solution y(t) = x(t + r). Notice that if k is 
small ('+' in Fig. 4(a)) it appears that (t) could even be larger than r. However 
this result does not take into account that the quality of the synchronization 
is poor in this case (it corresponds to the grey region near the vertical axis 
in Fig. 3) and that the standard deviation a is large (see Fig. 5) indicating a 
bad synchronization quality. Notice, finally, that for each value of k there is 
a maximum anticipation time, in agreement with the rather sharp transition 
between synchronized and desynchronized regimes shown in Fig. 3. 

Figure 4(b) plots the mean anticipation time (t) as a function of k for different 
values of r. The main result is that for each value of r there is an interval of 
values of k, K min < k < K max , such that (t) ~ r (the plateaus in Fig. 4(b)). 
For small values of the coupling, k < n m i n it is (t) ~ 0, and this reflects 
that the two neurons are synchronized (not anticipately) due to the common 
external forcing (this parameter region corresponds to the dark region close 
to the vertical axis of Fig. 3). If k, > K max the anticipation is lost due to the 
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chaotic behavior of the slave. 



a) b) 




Fig. 4. (a) Mean anticipation time as a function of the delay time r for the following 
values of the coupling strength: k = 0.15 (+), k = 0.25 (*) and k = 0.45 (O). (b) 
Mean anticipation time as a function of k for r = 1 (+), r = 2 (*) and r = 3 
(O). The results come from numerical integration of Eqs.(3-4) using the parameters 
indicated in Fig. 2. 
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Fig. 5. Plot of the standard deviation of the anticipation time, a, as a function of 
(a) the delay time r and (b) the coupling k in the same cases as in Fig. 4. 

More information about the quality of the anticipated synchronization is ob- 
tained by looking at the dispersion in the values of t™ — tj. In Figs. 5(a) and 
5(b) we plot the standard deviation, a, in the same cases as in Figs. 4(a) 
and 4(b). Of course, the best synchronization quality can be defined as the 
one with a small number of errors and a small dispersion in synchronization 
time. In this sense, one can see in Fig. 5(a) that a is an increasing function 
of r, indicating that the dispersion (and the quality of the synchronization) 
worsens for large r. Notice also in Fig. 5(b) that in the interval of coupling 
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strength where good synchronization occurs, K min < k < K max , a decreases 
significantly. 

In conclusion, we have studied numerically the regime of anticipated synchro- 
nization in coupled FitzHugh-Nagumo model neurons subjected to the same 
random external forcing. A difference with previous studies is that the an- 
ticipated synchronization manifold is not an exact solution of the dynamical 
equations. However, we have shown that the slave can predict the pulse fir- 
ing quite accurately. We have introduced the normalized number of errors 
and the mean anticipation time to measure the degree of synchronization. We 
have shown that the anticipation phenomenon is robust and exists on a wide 
parameter region. 

We would like to end by pointing out that synchronous neuronal oscillations 
underlie many cortical processes. Whether the results of this paper are of any 
interest to biological systems is open to speculation, but we hope that our 
numerical results will stimulate the search for anticipated synchronization in 
real biological neurons. 
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